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Abstract- In this paper we consider the two step process cubic convergent extrapolated Newton’s method 
which requires the evaluation of the first derivative and suggest a two-step process of the same method free 
from the derivatives by using backward difference approximation. It is shown that the new method has a cubic 
rate of convergence as that of Extrapolated Newton’s method and the efficiency index of this method is \pS . 
Few examples are given to illustrate the efficiency of this method compared to Newton’s method and 
extrapolated Newton’s method. 

Keywords'. Nonlinear equations ; Iterative method ; Convergence Criteria ; Newton ’s Method ; 
Cubic Convergence 

I. INTRODUCTION 

We consider finding the zeros of a nonlinear equation 

/O) = o (i.i) 

Where / \D ezR —>R is a scalar function on an open interval D and f (v) may be algebraic, 
transcendental or combined of both. 

The famous well known quadratic convergent Newton’s method for finding the root 'oc' of (1.1) is given by 

/OJ 

X n + l= X n~J^ (n = 0,l,2,„.) (1.2) 

Starting from the initial approximation X Q which is in the vicinity of the exact root, V.B. Kumar, Vatti et. al. 
[15] suggested a cubic convergent method by extrapolating the newton’s method (1.2) known as extrapolated 
Newton’s method which requires three functional evaluations i.e; f (v), /'(v) and /"(v) in each 
iteration, defined by 

- 

X n+1 = X n~ (n= °’ !’ 2 "') 0 - 3 ) 

Here 'ct' is the relaxation parameter and its optimal choice given by 

2 

a= (1.4) 

2 ~Pn 
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Pn 



jwtBmk 



(1.5) 



(n= 0, 1,2...) 



( 1 . 6 ) 



A variant of Newton’s method suggested by Fernando and Weerakon [13] which is a two step iterative process, 
defined by 

X n + l= X n-^—^—r~ 

f^n) + f i x n+ 1) 

„ * 

where X - = X„ 

n + f n fXx j 

This method (1.6) has a third order convergence. 

It is to note that the efficiency indexes of the methods (1.2), (1.3), and (1.6) are and respectively. 

In this paper, we present the method (1.3) as a two-step iterative process free from all the derivatives in section 
2. In section 3, the convergence criteria of the new method is discussed where as in the concluding section 
several numerical examples are considered to exhibit the efficiency of the developed method. 

II. DERIVATIVE FREE EXTRAPOLATED NEWTON’S METHOD 

The cubic convergent extrapolated Newton’s method free from second derivative suggested by V.B. 
Kumar, Vatti et. al. [16] is given as: 

For a given Xq , compute the approximate solution X n+ ^ by iterative scheme: 

fr*„) 1 



y n = X n ~ 



. . (n=0, 1,2...) 

[f^ n )-f(y n )\f\x n ) 

f(*n) 

f '(x n ) 

As it is known that the backward difference approximation for the first derivative for /'(x) at x is 

/(x)-/(x-A) 

nx)~ h 

Replacing 7? ’ by /(x^ ) , the backward difference approximation for the derivative / (x^ ) at X n is 

xj) 

/ (^n) 

Now the method (2. 1) takes the form 



( 2 . 1 ) 



n* n > 



X„, 1 = X 



f\x n ) 






where 



‘” +1 f(.x n )-f(x n -f(xj)lf(x n )-f(y n ) 

■A*„l 



(n=0, 1,2...) 



y n =X n~ 



We now define the following algorithm. 

Algorithm 2.1 : For a given Xq , compute the approximate solution X fl+1 by iterative scheme: 




V.B.K. Vattilet al. /International Journal of Applied Sciences and Innovation 



X n ,1 = X — 



f\ X n) 


f( < X n) 


f(x n )-f(x n -f(x n )) 





y n = X n~ 






f ( X n )- f ( X n~ f ( X nV 
This algorithm can be called as a two step derivative free extrapolated Newton’s method and requires three 

functional evaluations. The efficiency index of this method is >/3 . 



III. CONVERGENCE CRITERIA 



Theorem 3.1: Let (X G D be a simple zero of the function f D <Z R — » R for an open interval D. If Xq 
is in the vicinity of exact root (X , then the algorithm (2.1) has third order convergence and it satisfies the error 



equation. = 



--1 



e„ where C = X 



Proof: If (X be the root of f(x~)= 0 and 
e n =X n ~ a 

be the error at n' h iteration. Then by Taylor’s series, we have 

f(x n ) = f(a + e n ) 

= m + f\a)e n+ el f ^ +e l f ^ + e' I 

= c l e „ + c 2 e l + C 3 e l + C A + c S e l + 0<e n> 
Where 

c fJ W 

' J \ 



7 4! 



5 rcay 






(3.2) 



— ( c i 36*2^1 + c*2 + 3c^c^ + 2c^C2 2^ +^3)^ 

+ (c 4 + C 2 (C2 +2c 1 c 3 )-5c 1 c 4 -5c 2 c 3 +6c^c 4 -4c^c 4 + c^c 4 +6c 1 c 2 c 3 -3 c 2 c 2 c 3 )e^ +0(e 5 n ) 

(3.4) 

Subtracting (3.2) from (3.4) we get 

f{x a )-n x n -f{x B y) 

= c\e n + (-c 2 c 2 + 3 c 2 c 1 )e 2 + (c 3 c 2 - 3 c 3 c\ - 2 c y c\ + 4c 3 c 4 + 2c 2 )e 2 

+(-c 2 (c 2 + 2c 1 c 3 ) + 5c 1 c 4 +5c 2 c 3 -^c\c 4 ^\c\c 4 -c\c 4 -6c x c 2 c 3 ^3c\c 2 c^)e\ +0(e \ ) 

(3.5) 

Also [/(aJ] 2 = [< 7^ + c 2 e 2 + c 3 e 2 + c 4 e^ + c 5 e^ +0(e^)] 2 

= c i e fl + ^ c i c 2 e n (^2 + ^^1^3 )^/; +(2c 2 C 3 +2c 1 C 4 )e /J (3.6) 

From (3.5) and (3.6) we get 
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[/ Q „)] 2 

ftxJ-ft.X'-ft.xJ) 

q 2 e 2 + 2c 1 c 2 e i n + (c 2 + 2 c 1 c 3 )e^ +(2c 2 c 3 + 2 c 1 c 4 )e 3 +0(e 
{ c\e n + (-c 2 c 2 + 3c 2 c 1 )e 2 + (c 3 c 3 -3 c 3 c\ -2c 1 c 2 +4c 3 c 1 + 2c 2 )e 3 
+(-c 2 (c 2 +2c 1 c 3 ) + 5c 1 c 4 +5c 2 c 3 -6c 2 c 4 +4c 2 c 4 -c\c 4 - 6 qc 2 c 3 +3c 2 c 2 c 3 )e* +0(e^)} 



? °2 2 ( c 2 + 2 ^^) 3 ( 2 c 2 c 3 - 2 ^^) 4 



(-c 2 c 2 + 3 c 2 c 1 ) ( c 3 c 4 - 3 c 3 c\ - 2 c 4 c 2 + 4 c 3 c x + 2 c \ ) 2 



, + 2c i c s) + 5c 1 C 4 + 5c 2 c 3 - 6c i C 4 + 4c i C 4 - C 1 C 4 - 6c i C 2 C 3 + 3c i C 2 C 3 ) ]g 3 +0(e 4 } 



2 ( c 2 + 2 c i c 3 ) 3 ( 2 c 2 Cg + 2 c^c^) 4 

+ 2 — + — + — — 



± I (~ C 2 C 1 + 3 c 2 c l) c , ( C 3 C 1 -3c 3 q 2 -2qc 2 + 4c 3 q + 2c 2 ) ^ 2 



, (~ c 2 c i +3c 2 q) c2 | 2 (~c 2 c i +3c 2 c 1 )(c 3 c 1 -3^ -2^ + 4c 3 q + 2c 2 ) ^ 3 



(-c 2 q 2 +3c 2 q) ^ 2 (c 3 q 3 -3c 3 q 2 -2qc 2 +4c 3 q +2c 2 ) ^ 3 



, ( c 2 c 1 +3c 2 q) +2 £2_ e 2 2 c 2 ( c 2 c 1 +3c 2 c 1 ) ^3 | (c 2 +2c t c 3 ) ^ 3 + ^ 4 

4 n ^ n ^ 2 n 2 n y n) 



c 7 ) ? (3c^cf -Cl Co -2c 3 ! c? -2CoCi + 2cZc^ + c?c?) 3 , 

e n + | c 7 - 2 | e 2 + v 3 1 LJ LJ 3_J LJ 2_iZ e 3+0( e 4) 
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y n = X n ~ 



y n =« + 



f(x n )-f(x n -f<ix n )) 



e + | C C 2 | e 2 , ( 3 c 3 C 1 -q^3- 2c l C 2 ~ 2 c 3 C 1 +2c 1 C 2 + C 2 C l\ 3 



|c2 (3c 3 q 4 -cfo -2 Ci V 2 2 -2c 3 q 3 +2q 2 c 2 2 +c 2 2 q 4 ) ^ 3 | 



From (3.8), we have 



/ x 7 (3c,c, 4 -cfc, -2^^ -2c,c\ 3 + 2c, 2 c? +c?c 1 4 ) o 4 

/o^ ) = ( c 2 - c 2 c i ) e l - — 0 — 1 2 - 3 1 — ^ — yyy e \ + oo 4 ) 

c i 

C i e n +C 2 e n + C 3 e l + C 4 e t + c $ e \ + ) 



From (3.2) and (3.9) we have 

fixJ-fiyA 



c.e n +c,c 1 e n + 



= 1 + M--C, e - 



(3c 3 q c x c 3 qc 2 2c 3 c x +2qc 2 ) ^ 2 



Combining equations (3.7) and (3.10), we get 

life r -ife# 1 



+ 1 c , _ c jl ] e z + &£. - c i 5f 3 zM/i z/fA iMd u!f!i e 3 + 0(e 4 ) 



*(? 



1+ — — c 9 e - 



(3c 3 c 4 -cfc 3 -q 3 c 2 2 -2c 3 q 3 +2 c\c\) ^ +Q ^ 



with (3. 1 1), (3. 1) and (2.5) one can have 



e n + l +a = e n +a ~ \ e n 




which yields 



+0(e*) 

(3.8) 

(3.9) 

+ C> On) 

(3.10) 

(3.11) 

(3.12) 
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It shows the method (2.5) has cubic convergence. 

IV. NUMERICAL EXAMPLES 

We consider few numerical examples considered by Fernando and Weerakoon [13] and by Grewal [1] and the 
method (2.5) is compared with the methods (1.2), (1.3), (1.6) and (2.1). The computational results are tabulated 
below and the results are correct up to an error less than 0.5x 10’ 20 . 

Table 1: Numerical Comparison 



Function 


*0 


number of iteration for each method 


ROOT 


ax) 




1.2 


1.3 


1.6 


2.1 


2.5 




(1)at 3 + 4x 2 -10 


1.8 


6 


6 


4 


3 


3 


1.365230013414097 


(2) cos(x)-x 


4 


41 


7 


8 


5 


4 


0.739085133215161 


(3) (x - 1) 3 - 1 


1.8 


6 


4 


4 


3 


3 


2 


(4) e x_1 +x-3 


0 


8 


5 


5 


4 


4 


1.442854401002388 


(5) x log 10 X — 1.2 


0.5 

2.5 


7 

5 


7 

5 


6 

4 


4 

3 


3 

2 


2.740646095973693 


(6)2jr-log 10 x-7 


3 

4 


5 

4 


4 

3 


4 

3 


2 

2 


2 

2 


3.789278248444742 


(7) xe x - cos AT 


0 


8 


5 


6 


5 


4 


0.517757363682458 


(8) e x sin x-1 


-0.2 


8 


5 


6 


5 


5 


0.588532743981861 


(9)e x -1.5 -tan -1 at 


-7 


7 


4 


6 


4 


4 


-14.101269772739964 



CONCLUSION 



With the number of iterations and the root of the respective equation tabulated for each of the methods, we 
conclude that the method (2.5) which does not require evaluation of the derivatives of the function, has the same or 
better rate of convergence compared to the methods considered in this paper. 
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